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We discuss a new ringdown frequency mode for vacuum perturbations of the Kerr black hole. We 
evolve initial data for the vacuum radial Teukolsky equation using a near horizon approximation, 
and find a frequency mode analogous to that found in a recent study of radiation generated by a 
plunging particle close to the Kerr horizon. We discuss our results in the context of that study. 
We also explore the utility of this mode by fitting a numerical waveform with a combination of the 
usual quasinormal modes (QNMs) and the new oscillation frequency. 
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I. INTRODUCTION 

Black holes are born when a massive star exhausts 
its nuclear burning processes, leading to a runaway col- 
lapse where gravity dominates over all other interactions. 
They can also be produced by the merger of binary sys- 
tems containing compact stellar remnants, such as neu- 
tron stars or smaller black holes. Stellar collapse and 
binary mergers resulting in black holes are astrophysical 
processes where it is expected that gravitational effects 
are strong, resulting in regions of high curvature. Obser- 
vations of such processes would provide a strong test of 
General Relativity. 

Gravitational wave astronomy will provide a power- 
ful tool for investigating astrophysical processes involv- 
ing highly curved regions of spacetime. In the absence of 
external fields and matter, black hole binary mergers are 
completely invisible in the electromagnetic spectrum, and 
no light can reach observers from the interior of a mas- 
sive star undergoing core collapse. In these situations 
gravitational waves are expected to carry away as much 
as a few percent of the total mass energy of the system, 
and can provide direct information about these otherwise 
unobservable events (see e.g [l|). 

In this study, we focus on the gravitational wave sig- 
nal produced in the final stages of the birth of a black 
hole, when the gravitational waves can be described using 
linear perturbation theory on the Kerr spacetime. Mea- 
surement of such waves could provide a key test of the 
so called "No-Hair Theorem" of General Relativity. The 
No-Hair Theorem is the statement that stationary black 
hole spacetimes are described completely using only a 
few parameters, namely mass, spin, and electric charge. 
This theorem has been proved for the case of Einstein- 
Maxwell black hole solutions, through the uniqueness 
theorem for the Kerr- Newman black holes [2|. Thus, 
when a black hole is born in a merger or stellar collapse, 
the resulting object must radiate away all of its multi- 
pole moments ^ > 2 over the course of a ringdown phase. 
This phase involves emission of gravitational waves in a 
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well known spectrum of exponentially decaying frequency 
modes, called the quasinormal modes (QNMs) l3|-p|, and 
also late time "tails" which have a power-law fall-off in 
time. Observed deviation from QNM oscillation in the 
ringdown phase would be indicative that the spacetime 
is not represented by perturbations on a Kerr spacetime, 
and so would be a violation of the No-Hair Theorem [9- 
i. 

In addition to this test, detailed study of QNM ring- 
down is a key component in detection of gravitational 
waves in the first place. Accurate theoretical and numer- 
ical gravitational waveforms are necessary for the suc- 
cess of the method of matched filtering, which will be 
used to extract the faint signal from the noise in these 
experiments. Matched filtering uses a gravitational wave 
template to filter the noise and determine if the wave 
is present. Accurate modeling of the ringdown phase is 
then necessary to build useful theoretical templates. 

In this study we focus on black hole mergers, which 
provide a cleaner system with definite numerical predic- 
tions, and for which the possibility of detection is higher. 
The recent strides in numerical relativity [lol | have al- 
lowed several groups to solve the problem of binaryinsgi- 
ral and merger completely for the first time (see Il2j | 
for recent reviews) . Such simulations have provided enor- 
mous insight into binary mergers, and indeed they can 
serve as a test bed for the theory of black hole perturba- 
tions, in addition to providing complete theoretical grav- 
itational wave templates. However, the computational 
expense of such simulations prohibits their use in gener- 
ating a large bank of templates for use in matched filter- 
ing. As such, a three stage, semianalytic approximation 
scheme has been developed to treat binary inspirals. This 
method has the advantage of reducing the computational 
expense for template generation. Also, analytic methods 
help to build intuition into the physical processes of the 
merger. 

The first stage is the long, quasistatic decay of the orbit 
of the binary, which is treated using the Post Newtonian 
approximation to General Relativity. The next phase is 
the rapid merger of the binary, requiring full numerical 
treatment (though various methods have been employed 
to approximate the entire merger, e.g. L1.3.J JJj)- Finally, 
once the two compact objects are surrounded by a com- 
mon horizon, the system can be approximated by the 
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evolution of perturbations of the final Kerr spacetime. 
The radiation generated in this phase is governed by the 
Teukolsky equation [l^ . The QNM frequencies are given 
by the allowed spectrum of the Teukolsky equation, when 
physically appropriate boundary conditions are imposed. 
The QNMs are located at the poles in the Green's func- 
tion of the radial Teukolsky equation, and are found using 
a variety of methods (see e.g. [5] for a recent review). 

Generally, it has been assumed that the QNMs make 
up the entire spectrum of oscillations during the ring- 
down phase after merger. Here we seek a new frequency, 
characterized by the properties of the Kerr horizon. We 
are inspired by a study by Mino and Brink (l6j . which 
investigated the radiation of a point particle falling into 
a Kerr black hole, using a near horizon expansion to find 
the radiation analytically. As the infalling particle ap- 
proaches the horizon, its trajectory in Boyer-Lindquist 
coordinates asymptotes to pure angular motion around 
the black hole with frequency fin, 



2Mr, 



= M + y/AP - a2 



(1) 



Here a is the spin parameter, M is the mass, and r+ is 
the Boyer-Lindquist radius of the outer horizon. From 
the viewpoint of an observer at infinity, the particle is 
frozen at the horizon, corotating with it and sourcing 
radiation at its rotation frequency. Calculations by Mino 
and Brink show that the radiation arrives at future null 
infinity with an exponential decay. 
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Here ^4 encodes the out-going radiation, as discussed 
fully in Section H, and m is the azimuthal quantum num- 
ber of the radiation. The decay rate gn is the surface 
gravity, given by 



gH 
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(3) 



The frequency here does not depend on details of 
the particle's energy and momentum, because the par- 
ticle's late-plunge trajectory is essentially universal in 
the Boyer-Lindquist coordinate system. This suggests 
that this oscillation mode, the "horizon mode," may be 
more general than this single case considered by Mino 
and Brink. If we take the naive point of view that the 
late stages of the merger can be approximated by gravi- 
tational perturbations falling onto a final black hole, then 
we have a situation where the infalling perturbations will 
source out-going waves like point particles. Though this 
viewpoint is crude, it does suggest a search for this new 
frequency mode in post-merger ringdowns. 

In this paper we will argue for the existence of a hori- 
zon mode (HM) with a frequency of mrtn and a decay 
constant which we find to be an integer multiple of gn- 
We find that the particular decay rate depends on our 
model for how the spacetime transitions from the non- 
linear merger into the regime of first order perturbations 
on the Kerr spacetime. 



This paper is organized as follows: Section HA pro- 
vides a simple argument for the presence of this fre- 
quency. In Section IIB we derive this mode through a 
direct construction, using a simple model for the transi- 
tion from merger to ringdown. In Section IIC, we explore 
the consequences of a different model for the transition. 
In Section HI, we reconcile our results with those of Mino 
and Brink. In order to test the utility of this new HM, 
in Section IV we use the HM in combination with the 
QNMs to fit a waveform generated by full numerical gen- 
eral relativity, and compare fits that include the HM to 
fits with the QNMs alone. 



II. THE NEAR HORIZON APPROXIMATION 

We first present a heuristic argument for the presence 
of a HM analogous to that of equation ([2]) in the solu- 
tions to the Teukolsky equation. We then derive the HM 
by evolving initial data for the Teukolsky equation in a 
near horizon approximation. Finally, we investigate the 
consequences of a different model for the transition of the 
spacetime into the regime of linear perturbation theory. 

A. Simple Argument for a Horizon Mode 

In the Boyer-Lindquist (BL) coordinate system, the 
components of the Weyl tensor which represent out-going 
perturbations of the Kerr spacetime are represented com- 
pactly by the Newman-Penrose curvature scalar. 

Note that we use a metric signature of ( h ++) 

and use the appropriate conventions of p7| for Newman- 
Penrose (NP) qua,ntities such as \l/4. We use the Kinner- 
sley null tetrad 
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where 



(5) 
(6) 

(7) 

(8) 
(9) 
(10) 

and the overbar represents complex conjugation. 

With these choices, ^4 satisfies a separable linear wave 
equation [l5j . and can be written as 

^^{t,r,eA)^ j dw^e-*'^*+™^i?,„<,(r)5,™^(0) . 
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Here, SimuAd) are the spin- weighted spheroidal harmon- 
ics [1^, [231 , with the appropriate spin weight for ^'4, 
s = —2. In the limit lu ^ 0, they reduce to the spin- 
weighted spherical harmonics ^2Yim{0,(j)). The radial 
function Rimuj(r) is the solution to the radial Teukolsky 
equation [1^. Note that ^4 vanishes in the background 
Kerr spacetime, and its perturbed value is independent 
of tetrad perturbations and gauge transformations. 

We can see the relationship between the scalar ^"4 and 
the out-going gravitational waveform via its asymptotic 
form near future null infinity. In this limit, for asymp- 
totically flat spacetimes (see, e.g. [2ll|). 



*4(r 



3) = -dtih+ -ihx) 



(12) 



Here the + and x indicate the polarization of the gravi- 
tational waves. 

In the BL coordinates, using the Kinnersley tetrad, the 
asymptotic behavior near the horizon of the two homo- 
geneous radial solutions are 



e''^'" out-going 



(13) 



m-gomg 



with k — Lo ~ mfln- Together with the separation of 5*4, 
equation (jlip. these solutions are associated with out- 
going and in-going radiation at the horizon. The tortoise 
coordinate r* is defined by dr* /dr = {r^ -f a^)/A. Note 
that r* — > —00 as r — >■ r+. 

One of these two solutions is selected out as unphysi- 
cal, based on its behavior near the event horizon. Here 
we repeat an argument first presented by Teukolsky [l^ . 
We demand that fields neither vanish exactly nor diverge 
at the horizon when measured by a physical observer. 
Near the horizon, the trajectory of any freely falling ob- 
server will approxima tely match that of an in-going null 
geodesic (see e.g. [H,|13): independent of the observer's 
energy or angular momentum, with 



t 



= const, 6 = const, </> = flHt . 



(14) 



This means that for distant observers, the infalling tra- 
jectory does not appear to enter the horizon, but instead 
asymptotes to it while circulating around the horizon 
with a constant frequency. A more natural set of co- 
ordinates is in-going Kerr coordinates, {v,r,9,(j)), where 
dv = dt + dr* , and dip = d(j) + a dr* /{r'^ + a^). In-going 
null geodesies lie on lines of constant v and (p. 

In this in-going coordinate system, the metric does not 
become singular at the horizon. However, even in these 
coordinates, the Kinnersley tetrad used to define ^4 be- 
comes singular at the horizon. This can be repaired by 
using a Lorentz transform to boost into the reference 
frame of an infalling observer who carries a nonsingular 
tetrad, namely 
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Here the subscript "in" indicates the regular in-going 
tetrad. The tetrad, written in the in-going Kerr com- 
ponents, is 
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1, 



::0, 
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With this tetrad, the physical observer measures a cur- 
vature scalar of 



2(r2 + a2) 



n 2 



A 
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(17) 



As a consequence, the two radial solutions in equa- 
tion (HH) correspond to *f - ^-tu;v^t7n4> ^in _ 
^-2g-iw'Ug4m0g2ifcr*_ rpj^^ sccoud diverges at the hori- 
zon, and so is selected as unphysical. In other words, 
if waves emerge from the horizon, then the in-going ob- 
servers will see a diverging curvature due to blueshift 
effects. Note that while it is a particular observer that 
carries the tetrad of equation near the horizon, these 
general results hold for all physical observers. This is be- 
cause the tetrad that another physical observer carries 
can be related to that in equation (jl6p through nonsin- 
gular Lorentz transformations. 

However, if the frequency lo is complex, this divergence 
can be removed. To see this, we let lo = mfln — 17, and 
seek an appropriate 7 (the real part of the frequency is 
chosen so that the observer does not measure increasingly 
rapid oscillations when approaching the horizon). We 
note that near the horizon, A ^ e'^^'^^ . For the physical 
observer, then, 



exp —iuv + imcj) — Agnr* + 2'yr* 



We see from this that if 



7 = 2g// , 



(18) 



(19) 



this solution remains regular at the horizon, and the so- 
lution decays along the observer's worldline in just such 
a way that the growth is compensated for. This par- 
ticular frequency is selected out by a physically allowed 
solution of the Teukolsky equation, and so we must con- 
sider its place in the usual QNM spectrum. Modes with 
7 > 2gH decay exponentially for the physical observer. 
In this sense, the decay rate of equation gives the 
least damped, physically reasonable mode in this simple 
argument. 
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B. The Near Horizon Limit: 
Model 



The E Boundary 



In order to evaluate this expression, we insert the ex- 
phcit form of the radial Green's function, 



With our heuristic argument in hand, we now derive 
the HM using a specific ansatz. In order to keep our 
results relevant to the problem of compact binary inspi- 
ral, we consider the following model. We imagine that to 
the future of some spacelike hypersurface S, the space- 
time can be described with linear perturbation theory on 
the Kerr spacetime, while to the past of S the space- 
time may be nonlinear. We denote E by setting the 
Boyer-Lindquist time coordinate t = 0. The past of S 
represents the inspiral and merger phases of binary co- 
alescence. To the future of S we can use the Teukolsky 
equation to evolve initial perturbations on S forward. A 
similar ansatz has been used in the Close Limit approach 
[23{ , and Lazarus project [13] which used numerical inte- 
gration of the Teukolsky equation to evolve initial data 
on an initial time slice [2J, [2^ . We refer to our model as 
the E boundary model. 

Given this ansatz, we evolve initial data V'|t=o = 
/o~'*^'4|t=o and dtip\t=o, using the Green's function. The 
full details of the analysis are presented in Appendix A. 
A near horizon expansion allows us to obtain part of the 
evolution analytically. Physically, we postulate that just 
after the merger of the binary, the perturbations are con- 
centrated in a small region near the horizon, so that the 
initial data used in the Green's function evolution only 
has support in a small region near the horizon. This 
expansion allows us to use the asymptotic form for the 
Green's function, and to keep terms only to first order 
in e = (r — r_|_)/r_|_ ^ 1. While we focus on the physical 
picture where the perturbations are concentrated near 
the horizon, our results hold for the evolution of the ini- 
tial data which is near the horizon even if the data on 
the initial surface extends to large r. In addition, due 
to redshift effects, this finite region near the horizon pro- 
duces (decaying) radiation over an infinite region of null 
infinity. 

Specifically, let ip\t=o and dtip\t=o be nonzero only be- 
tween r+ and (1 -I- ^)r+, with ^ ^ 1. We truncate the 
integrals of the Green's function over the initial perturba- 
tion to this small region. To first order in distance from 
the horizon, A 2AIr^Ke, with k = \/l — a? /A-P. From 
Appendix A, equation (|A23p . we have that \['4 takes the 
form of equation (ITTT) with the radial function Rimi^{r) 
given to leading order in e by equation (IA24[) . 



dr' 



2Mr+(KtY 
{2Mk + ima)aemu]{r') 



2(Mr+)2(Ke)3 



Glrauj{r, r') 



(20) 



where the functions aimujiT) and l3imuj{r) can be found 
from the initial data, using equations (jAlsp . (|A16|) . 
(IMT|) . and (|M2l) . The function G£™„(r,r') is the fre- 
quency domain radial Green's function. 




< r 
> r 



(21) 



The functions i?"^^^ and Rf^^^^ are two homogeneous solu- 
tions to the radial Teukolsky equation, with the up-mode 
(no radiation from past null infinity) and in-mode (no 
radiation from the past horizon) boundary conditions, 
respectively. They have the asymptotic forms 



^ — ikr* 



(22) 



R 



up 



ikr* 



/^trans„3 , 



(23) 



where the Wronskian Wemoj is given by Wtn 
2iujB\ 



Since we are interested in the waves at infinity, r — > 
c», we insert the appropriate asymptotic expression for 
Also, since the integral extends only over the 



R 



near horizon region, we insert the asymptotic expression 
Thus, 



Dtrans A 2 „ 



ptrans 



2iujB 



(24) 



imuj 



dr € ^e' 



ikr* 



-{iMujr^ + 2Mk + ima)aimLo{r) 



(25) 



In order to complete the integration, we need to know 
how the initial data behaves to leading order in e. We can 
write to leading order aimujir) ~ and (3imuj{r) « 

(^irnuj^'^- We find the leading order powers n and p by 
essentially repeating the argument given in Section IIA, 
with some additional care. 

Here we cannot follow a single observer who falls past 
the horizon, since we wish to know the behavior of ^4 
on the initial slice S. We consider instead a family of 
accelerated observers who cross the surface t = at a 
variety of radii, extending all the way to the horizon but 
not penetrating it. We choose the four- velocities of these 
observers at this initial surface to vary with r, u(r), such 
that all members of this family carry the regular tetrad 
P6I) . These observers measure an initial perturbation 
in their own frames, and for an observer at some r the 
measured perturbation is related to that expressed in the 
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Kinnersley tetrad by = [2{r'^ + a?)//^]'^'^^ ^ e"2*4. 
Since all the observers are physical, we expect that the 
perturbation they measure does not diverge as we take 
the limit r — >■ r+, moving along the family of observers. 



However, this requires that ^4 ~ 
n = 2, p = 2 on this interval. 
Inserting this into we have 



e . Thus, we require 



Mr +13^^^ + {iMujr+ + 2Mk + ima)a°^^^ 



'■+(1+?) 



dr 



-ikr* 



Mr+I3'^„^^ + {iMujr+ + 2Mk + inia)a°e„ 
2 



(26) 



Focusing on the integral, which we denote /, we have 



cj - {mSXH - 4:igH) 



j2-ik/2gH)\n(, 



(27) 



Here, there is a pole in the denominator which will select 
out the frequency 



(E Model) 



(28) 



when equation is integrated over uj. 

Combining equations (HH), ([Ml), dM]), and (|27l) . we in- 
tegrate over oj. We close the contour in the lower half 
plane, selecting out the poles of equation ([M)) by the 
residue theorem. The zeros of -B"^^ comprise one set of 
poles in the lower half plane, and these poles give the 
usual QNM frequencies. We wish to focus on the con- 
tribution from the additional pole at loh, and so from 
here we ignore the poles coming from B^^^. Similarly, 
we will not consider here the influence of the pole at 
a; = (actually part of a branch cut along the negative 
imaginary axis, which generates the late-time power-law 
tails). Also, we will not consider the possible poles in 
the terms 01^^^ and which receive their frequency 

dependence from the projection of the initial data onto 
the spheroidal harmonics. 

This integral converges for t — r* + hi^/2gH > 0; 
otherwise, we must close the contour in the upper half 
plane and the integral vanishes. Noting that as r — >■ 00, 



P 



we have 



*4 



\ ^ 7 —iuJH{t — r*)+imcji q / a\ 



hn 

xh( t-r* 



In^ 



(29) 



{Mloht 



ma — 2iMK,)a 



,0 



(30) 



where H is the unit step function. In the above, all fre- 
quencies are to be evaluated at ojh from (^5)) . As <^1, 
the waves at infinity appear at late retarded times. This 
sharp turn-on of the wave is an artifact of our trunca- 
tion of the integral at (1 -|- ^)r+. A smoother falloff of 
the initial data with increasing radii would result in a 
smoother turn on of the wave at infinity. These waves 
at early times are sourced by initial data on E which 
cannot be evolved using the near horizon approximation. 
We see also that the waves continue to reach infinity for 
all retarded times after the turn on. As mentioned pre- 
viously, this is due to redshift effects near the horizon, 
which stretch the radiation from the finite near horizon 
region out over an infinite region of null infinity. 



C. The Characteristic Boundary Model 

While the frequency of the radiation in equation ^7} 
matches the result of our heuristic argument, the decay 
rate does not. The decay rate is determined by the radial 
behavior of the initial data, as we can see in equations 
(1^51) and (l?7l) . In order that luh = rnilH — '^Wh, the 
initial data would need to behave as ^"414=0 ~ which 
we have argued against based on our requirement that 
physical observers near the horizon measure nonsingular 
initial data. However, a change of our ansatz shows that 
the initial data can be proportional to e and still represent 
physical perturbations. In this case the surface bounding 
the regime of linear perturbation theory is an in-going 
null surface, instead of a surface of constant t. This differs 
from the E boundary model, and so differs from the Close 
Limit Approach. We will refer to this second model as 
the characteristic boundary model. 

As the spacetime transitions into the linear regime, the 
nonlinear perturbations radiate away towards infinity or 
down into the black hole. We imagine that the regions 
of nonlinear evolution are bounded by characteristics of 
the linear wave equation. This is a more physically mo- 
tivated assumption than a transition in spacetime prop- 
erties along the entire surface E. The in-going character- 
istics have a trajectory r{t) = r+(l -|- e"^^"^*"*"^) [l6|, 
where t — to ^ 1. The comparison of these two models 
in both BL {t,r) coordinates and the tortoise {t,r*) co- 
ordinates is given in Figure [1] In the {t,r*) coordinates, 
we see that the horizon is pushed to r* — > —00., and that 
the initial data of the E model is stretched out onto an 
infinite interval in r* . We see in both figures that the 
horizon is hidden behind the boundary characteristics, 
and our previous argument for the radial dependence of 
the initial data on E no longer holds. We must find a 
new way to determine the r dependence of initial data in 
this model. 

We will again evolve initial data on the constant time 
slice E, this time with support only on a small in- 
terval outside the boundary characteristics. At time 
t = 0, we set the inner boundary of the initial data 
set to be at r+(l -f ^1), and the outer boundary at 
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FIG. 1: (color online) Comparison of the E and characteristic 
models. In both figures, the trajectory of the observer dis- 
cussed and the characteristic boundary surface are plotted, 
with (9, (j)) suppressed. The thick line is the interval where 
the initial data is nonzero, and the two points are and ^2, 
respectively. The shaded regions correspond to the nonlinear 
regime of each model. Top: Comparison in BL coordinates 
of the two models, with the event horizon illustrated. Bot- 
tom: Comparison in {t,r*) coordinates, where the horizon is 
pushed to r* — >■ oo. 



r-|_(l + ^2)- We will use the same physical observers as 
before. To first order their trajectories are lines of con- 
stant V and 0, just like the in-going boundary character- 
istics. The physical observer who passes r-|_(l + ^2) at 
t = has a trajectory r'obs(i) = + e-29«(*-*ob=)). 
At i = the observer measures an initial perturbation 
^'Jj" ^ f^^^4(robs(0)). At a later time t, the observer 
measures ^4" ~ ^^^e'*^"*^4(robs(0)i and so the mea- 
sured perturbation grows exponentially in time. How- 
ever, the perturbation will also decay in time due to its 
evolution. We insist then that the decay be such that this 



observer (and similarly, all of the physical observers near 
the horizon) do not experience an exponentially growing 
perturbation. This will set the behavior of ^'4 on the 
initial surface. 

We examine then the perturbation as measured by the 
observer along his trajectory. Given data that behaves 



and Pi 



/39 e" on the initial 



surface, we again combine equations (pn)) - this time 
taking the asymptotic limits as r' — > r+ and r ^ We 
focus on the out-going solution only, since the observer 
will not measure the in-going waves. We find that 



uj — {mVLH — inign) 



X I e '23ff (t-tobs)(«-ifc''+/K) _ gln^i(n-ifcr+/K 



(31) 



Integrating equation with this radial function, and 
noting again that for the observer Vahs = t -\- r* is con- 
stant, so that t = Vohs ~ f*, we have 



*4(f'obs) ~ e 



--imQ,Ht—^ngHt 



H{t+t 



obs 



-H 2t 



In a 

2gH 



(32) 

For n — I, the decay of the perturbation along the world- 
line of the observer is just enough to cancel the exponen- 
tial growth. The initial data on E can be taken to be 
proportional to e. 

With the initial data, we can return our attention to 
the perturbations measured at infinity. Repeating the 
analysis of Section IIB with this initial data, we have for 
r ^ 00 



*4 



E 



H t-r* 



In 6 



-H t-r* 



In 6 



(33) 



with Zimuj as in equation pop but with uh given by 

ujh ~ miln — 2igH ■ (Boundary Model) (34) 

We have recovered the decay rate indicated by the heuris- 
tic argument of Section IIA. The difference of step func- 
tions here and in equation ([32]) is again due to the sharp 
truncation of the integral at each end of the interval on 
S. Now the radiation turns off due to the truncation of 
the initial data at the boundary characteristic ^ . At this 
retarded time the gravitational radiation would give way 
to radiation sourced by the perturbations in the nonlin- 
ear region of spacetime on S. For this second model it 
seems that a method for evolving data along character- 
istics would be better suited than evolution from a con- 
stant time slice. Such a characteristic evolution has been 
presented for the Schwarzschild black hole, for example 
in 26]. Others [27*1 have presented numerical evolution 



7 



of characteristics, again for the Schwarzschild black hole. 
Another possible formulation which would be natural in 
this context would be the use of an asymptotically hy- 
perboloidal spacelike surface in place of E, as discussed 
in m. 



III. RECONCILIATION WITH THE 
MINO-BRINK MODE 



In Section IIB and IIC, we saw that the condition 
that physical observers measure regular curvature near 
the horizon determines the decay rate for the gravita- 
tional radiation at infinity. In fact, the mode of Mino 
and Brink, with frequency w — ttiVLh — ign, has a decay 
rate which violates the regularity conditions discussed in 
both sections. Its decay rate is too small, for example 
compared to equations (p8)) and (|34)) with d ecay rates 
7 > 2(7i/. Thus, although the mode found in [1^ moti- 
vates our study, the two results are in disagreement. In 
this section, we first provide a simple alternative estimate 
for the expected decay rate of radiation from a point par- 
ticle, using the notation of Newman and Penrose for con- 
venience [29j . With this estimate as motivation, we then 
find that correction of an error in unexpectedly leads 
to the vanishing of the first order mode discussed there. 
We conclude that the actual leading order radiation from 
an infalling point particle has a decay rate which matches 
our characteristic boundary model in Section IIC. 



Point Particle Radiation in Newman-Penrose 
Formalism 



We wish to calculate the radiation generated by a 
point particle at the last stages of its plunge into a Kerr 
black hole. Once again, the radiation is described by 
the Teukolsky formalism, this time with the appropriate 
source term for the matter content. However, it will turn 
out to be convenient to make the near horizon expansion 
in the Newman-Penrose formalism, in order to obtain an 
estimate for the behavior of the radiation. As we show 
in Appendix B, we can write to leading order in A (e) 
near the horizon, using the Kinnersley tetrad. 



where 



In addition, we can approximate 



+ 47£))«'4 « 47rT4 , 


(35) 


D = l^V^ , 


(36) 


A = n^V^ , 


(37) 


r - M 


(38) 


^ ~ 2^ ' 



D 



A 



2A 



21] 



I 

-{dt - dr' + ^Hd^) = 



2A 



2S 



-i+,(39) 

J 

-L_,(40) 



and write 

M'^r\{L+L_ + AgHL+)^i « 47rI]AT4 . (41) 

In absence of T^^ this directly gives the asymptotic in- 
going and out-going waves of equation Let us spe- 
cialize 4'4 to a particular azimuthal quantum number, to, 
and we have 



L± — dt i: dr, + imU.H 



(42) 



We turn now to T4. We expand the stress-energy ten- 
sor of a point particle in terms of the azimuthal quantum 
number to to match the expansion of ^'4 implicit in (j42l) , 
which gives 



1 



9 

uf^u" dr* 
M*S dr 



S{r - r{t))S{9 - 6l(t))e™('^-'^(*» (43) 



■S{r* - r*{t))d{x ~ x(0)e""''^"'^'*^^ 



(44) 



Here we have defined x = cos 9, and used the properties 
of the delta function. The mass of the particle is given by 
/Lt. If we define v — t + r* and u = t — r*, and once again 
use the properties of the delta function, we can write 



S{r* - r*{t)) = S{ — : r 



du/dt 



(45) 
(46) 



where VQ{t) is the value at which the argument of the 
delta function in (|45p vanishes, which is at first order 
simply the value uq = const to which the trajectory 
asymptotes. Also, to leading order the trajectory will 
have t — —r*, so du/dt = 2 to leading order. 
We have then that 

(47) 

We must now project T^'' onto the null Kinnersley ba- 
sis in order to find T^, see (IB2p . This will result in a 
projection of the four velocity onto the Kinnersley basis, 
for example with u„ = Uf^n'^, and in this basis some of 
the components are vanishing or divergent as r — > r+ (as 
seen in Section HA). In order to examine the leading or- 
der behavior near the horizon, it is then best to express 
the four velocity components in terms of the regular, in- 
going basis, related to the Kinnersley tetrad by (ITSI) . In 
the in-going basis, we have 



and where the source term T4 is given by (IB2[) . 



(48) 
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with ttjj^"'"''™ all smooth and finite throughout the trajec- 
tory — including on and within the future horizon. The 
four velocity in the Kinnersley basis can be expressed as 



u = m'Ik + u"nK + u"mK + u'"mK 



(49) 



Ik + ^ riK + u[nmK 



2(r2 + a2) 



A 



(50) 



Finally, we lower the tetrad indices on the components 
of the four velocity using the null metric 



Vab = 



0-100 
-10 
1 
10 



This gives us 



ii 

m\2 



(51) 



(52) 
(53) 
(54) 



Let us first consider Tmrn, which gives the dominant 
contribution T^mm to T4, 



-(A + 27)Ar^ 



(55) 



It turns out that to first order ATmm = 0. We can see 
this by noting that at first order 



Cj.2 I „2^|2 



(56) 



(n,m\2 im{tti-nHt) 

- ^A;;;(i.-^^.L.) ^(--")^(^-^°)(^^^ 

with xo = cos 60 the value of x to which the particle 
asymptotes to at the horizon. Acting on this with A 
from (|40|) . we get 



AT- - + 



2S 



So we see that 



(58) 



(59) 



Physically, this is because A takes its derivative almost 
along the direction of motion, along which does not 
change to first order. This means that the exact con- 
tribution of Tfam to T4 must bc reexamined with other 
terms included, and its contribution is in fact at the order 
of Tffin- Thus, we expect 



(L+L_ + 4.gHi+)*4 ^ A2(5(w - uq) 



(60) 



We have a simple scenario: if we integrate across the 
V = vq surface, removing the derivatives dy = {dt+dr-) /I 



from the left-hand side, we will have a u-dependent ^'4, 
which obeys 



A^ 



-ignt 



(61) 



Here we have recovered the decay rate near the horizon 
discussed in Section IIC. This indicates that the corre- 
sponding decay rate of the waves as r — >■ 00 is that of 
equations However, this argument lacks the 

detailed calculations of Mino and Brink, who found a 
mode with a slower decay rate. Under examination, how- 
ever, it is an error in [l6[ which leads to a mode with a 
spuriously low decay rate. We discuss this in the next 
section. 



B. Eliminating the Leading Order Frequency Mode 
of Mino and Brink 



We turn now to the study by Mino and Brink, which 
we abbreviate as MB. In this study, the source term is 
evolved using the Green's function method much as Sec- 
tion II of our study, resulting in an integrand for the 
integral (I2ip which has a pole at w = m^ln — igu, i-e. 
equation (MB 3.7). The pole then selects out this oscil- 
lation frequency for the out-going radiation at infinity. 
However, it turns out that this pole is canceled out by 
terms in the amplitude Zi„iuj , when an error in MB is cor- 
rected for. We find this error in going from (MB A14) to 
(MB B4) . The first equation gives a piece of the Fourier 
decomposition of the source term T4, and is drawn from 
[30| . equations (2.21) and (2.25) therein. 

Taking the leading order contribution of (MB A14), we 
find that (MB B4) should read at leading order 



A__ 



mmO 







K 





,r+-iacosOo -200 -1 
r+ + tacosuo 



(62) 



where 60 is the value of 6 that the point particle asymp- 
totes to on the horizon, and -Bisco contains information 
on the particle's constants of motion. This equation dif- 
fers from MB by the factor of —1 in front of the term 
ikrjf. / K. When this difference is accounted for, we have 
for (MB 3.6) 

Zimu, cx (w - + igH){t^ - miln + ^ign) ■ (63) 

The first root of Zim^ then cancels the pole in the de- 
nominator of (MB 3.7). It also appears that the second 
root removes the frequency mode our rough argument in 
Section III A suggests. However, (|63p holds only to first 
order, and at second order there are additional terms not 
proportional to these roots. While a second surprising 
cancellation can only be ruled out by a careful study of 
the MB analysis at second order (or an equivalent for- 
mulation), it would seem unlikely that the next order 
of frequency mode would vanish as well. Such a careful 
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FIG. 2: (color online) A comparison of the real and imag- 
inary parts of cjh (dashed) and the first four QNMs, with 
^ = 2,m = ±2 (solid). 
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FIG. 3: (color online) The section of the numerical waveform 
given in Scheel et. al. over which to ranges in the overlap 
integral. Here we give both the real (solid) and imaginary 
(dashed) portions of Mor'i'i. 



study goes beyond the scope of this paper. However, in- 
vestigation does indeed show that at second order there 
is a HM with a decay rate of 2gH [3l| . 



IV. NUMERICAL STUDY 

In the previous sections we argued for the presence of 
a HM in the ringdown spectrum. In this section we test 
a numerical waveform for evidence of this mode. For this 
study we use the publicly available waveform generated 
by Scheel et. al. [32| by the evolution of an equal mass 
black hole binary through merger and ringdown. First 
we compute the overlap between the final portion of the 
numerical waveform and either a combination of QNM 
oscillations, or a combination of QNM oscillations and 
the HM, in order to see if the given combination is a 
good fit to the waveform. Next, we use a best-fit of the 
overlap to extract the mass and spin of the black hole 
from the waveform. Again, we compare the extraction 
using just the QNMs with an extraction which includes 
the HM. 

For our HM, we focus on the less damped of the two fre- 
quencies discussed in Section II, loh = mO.H — '^igH- This 
mode agrees with the frequency mode from the point par- 
ticle plunge discussed in Section HI, and also from the 
more physically motived model of the perturbed black 
hole's transition into the linear perturbation regime, dis- 
cussed in Section IIC. Throughout this section we will 
discuss only the dominant £ — 2, m — 2 mode in the 
spherical harmonic decomposition of the waveform. Note 
that, because of the rotation of the Kerr black hole, there 
is a Zeeman-like splitting of the QNMs into modes which 
corotate with the black hole and modes which counter- 
rotate with the hole. The counterrotating mode frequen- 
cies are equivalent to the £ = 2, m ~ —2 QNM fre- 
quencies, but with a negative real part of the frequency, 
^nhn — ~^ni-m (see iTj] for a detailed discussion). Each 
additional overtone we consider in the numerical analysis 
in the following sections thus adds two distinct modes. 

Figure [2] compares the real and imaginary parts of lOu 
with those of the first four QNM for € = 2, m = 2 
(corotating modes) and € = 2, m = — 2 (counterrotat- 
ing modes), as a function of g/M. These QNM fre- 
quency values are drawn from j33| and calculated us- 
ing the methods discussed in f5l|, whose values are used 
throughout this study. Note also that while many stud- 
ies refer to n = as the slowest decaying QNM, here we 
count overtones from n = 1. The corotating QNMs gen- 
erally have a higher frequency than the HM, and the first 
two QNMs have slower decay rates. The third corotating 
QNM has a comparable decay rate, and the fourth de- 
cays faster than the HM. Meanwhile, the counterrotating 
QNMs decrease in frequency with increasing a/M , until 
they oscillate slower than the HM. The decay rates of the 
counterrotating QNMs also remain nearly constant over 
the whole range of a/M, and so in this case the third 
and fourth QNMs both decay more quickly than the HM 
over a large range of a/M. As a/M goes to zero, luh 
ceases to drive oscillations, and the HM simply decays 
exponentially. Also, we see that as a/M goes to zero, 
the degeneracy between co and counterrotating modes is 
restored. Finally, in Table IIVI we give numerical values 
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(n, I, m) 
(1,2,2) 
(2,2,2) 
(3,2,2) 
(4,2,2) 
(5,2,2) 
(1,2,-2) 
(2,2,-2) 
(3,2,-2) 
(4,2,-2) 
(5,2,-2) 



MtJQNM 

0.52670 + 0.08129j 
0.51486 + i0.24581i 
0.49296 + 0.41513i 
0.46387 + 0.58873i 
0.43291 + 0.76035i 
0.31072 + 0.08874i 
0.27312 + 0.27733j 
0.21198 + 0.49963i 
0.15865 -H0.75811i 
0.12707 4- 1.030311 



MUJE 



2 0.37177 + 0.43089i 

TABLE I: Values of the QNMs and HM evaluated at the spin 
parameter a//M/ = 0.68646 appropriate for the final black 
hole of fsl. 



for the QNMs, evaluated at Uf/Mf for the final, merged 
black hole vi^hose spectrum we study [Hi (in units of the 
initial ADM mass of the binary), as well as the value of 
the HM at this spin parameter. 



A. Comparison of Overlaps 

Consider the final portion of the numerical waveform, 
tp{t)H{t — to), with H{t) the unit step function, and to 
a constant which we consider to be the time where the 
ringdown phase begins. We wish to see how well a wave- 
form made from a linear combination of damped sinu- 
soids, Vfc = e(~''^'=~^'=)(*~*'')i?(i - to), can be made to fit 
■0. Since we do not know a priori at what point in the 
numerical waveform the ringdown phase begins, we vary 
to as a free parameter in our study. This allows us to see 
where in the waveform our combination of sinusoids fails 
to be a good fit; at sufficiently early to we do not expect 
a particular combination of to model the chosen sec- 
tion of the waveform accurately. However, a combination 
of V'/c that fits the waveform well over a range of to that 
includes the early parts of the ringdown more accurately 
represents the frequency spectrum of the ringdown than 
another set of damped sinusoids that first fails at a larger 
value of to. 

For two waveforms si(t) and S2(t) cut off at to, we first 
define the inner product 



(Sl|s2) 




3920 3940 3960 3980 

to/Mo 



4000 4020 4040 



FIG. 4: (color online) Comparisons of 1 — Pmax for overlaps 
using just the QNMs (solid) and those which replace one QNM 
pair with the HM (dashed) . The top panel compares the first 
QNM pair (both co and counterrotating) with the HM alone. 
The second panel compares the first two QNM pairs with the 
first QNM pair and the HM. Subsequent panels compare the 
first n (n = 3, 4, 5 )QNM pairs with the first n — 1 pairs and 
the HM. 



the combination ^ akipk 



P 



2 _ k _ 



(65) 



where = {tpk\'<P) , -B/„i = (i^iltpm)- This maximum 
overlap characterizes how well the ipk can be made to 
approximate ip, given the optimum choice of ak- The 
maximization yields^ 



Pmax[V';{Wfc,7fc}] = 



(010) 



(66) 



Using this equation to compute Pmax, we take the 
first n QNMs (recall that each overtone includes two fre- 



si(t)s2(t)cit 



(64) 



The overlap, p, of two waveforms is given by the magni- 
tude of the normalized inner product of the waveforms. 
Our goal then is to maximize the overlap of 0(tr) and 



^ The Lagrange multiplier method yields AA^a = ABa, with A 
already the extremum. This means A should be the maximum 
eigenvalue of M = B~^y4At. However, since M only has one 
non-zero eigenvalue, we have A = trM = AtR^^A. 
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quency modes, both the co and counterrotating modes) 
and compute the maximum overlap with the numeri- 
cal waveform as a function of the starting time Iq. We 
then find the maximum overlap using the first (71 — 1) 
QNMs and the HM. Here, and throughout this section, 
we normalize our units by Mq, the sum of the initial 
ADM masses for the two black holes that merge [s^ . 
In these units, the mass Mf of the final black hole is 
given as Mj/Mq = 0.95162 ± 0.00002, and its spin a/ 
by af/Mo = 0.65325 ± 0.00004. For large values of to, 
only the least damped mode contributes to the waveform, 
and so any tpk that includes the first QNM will provide 
Pmax ~ 1- Therefore it is useful to investigate 1 — Pmax- 
Figure[3]plots the segment of the waveform over which we 
range to for the overlap calculations. Figure H] compares 
1 — Pmax on a log scale for overlaps using the QNMs alone 
to those including the HM. In the context of gravitational 
wave signals, it is preferable to have a good overlap for 
to close to the peak of the signal. This is when the grav- 
itational wave signal is strongest, and also the point in 
the waveform when the number of QNMs that make a 
significant contribution to the waveform is the greatest, 
before the most-damped QNMs become negligible. 



We see that the horizon mode alone provides a poor 
fit for the waveform for all values of to, and that the first 
two QNMs provide a better fit at earlier to than the first 
QNM and the horizon mode. Replacing the n"^ QNM 
with the HM gives comparable results at 71 = 4, and 
provides a (very slight) improvement at early to for n > 5. 
We note that this improvement becomes apparent at a to 
earlier than the peak value of \ip\, which occurs at tmax — 
3953. 8A/0. This means even for n > 5, the HM makes an 
improvement only when portions of the waveform which 
should not be well modeled by a set of damped sinusoids 
are included in the overlap. In fact, we find empirically 
that adding modes with low decay rates always tends to 
improve overlap calculations at early values of to. This 
is due to the fact that less damped modes will better 
fit the region near the peak of the waveform. In this 
case we would expect comparable results from the n = 3 
mode addition and the HM mode addition, and for the 
HM to improve the overlap compared to n = 4; however, 
the single HM must compete with the pair of co and 
counterrotating modes that make up the next overtone. 
For these reasons, we find that overlap comparisons do 
not provide a compelling case for the presence of the HM 
in the waveform, nor do they rule the mode out. 



As a second test, we would like to investigate the use of 
the HM in performing parameter extraction from a ring- 
down waveform. This practical test of the utility of the 
HM is more physically motivated than overlap compar- 
isons, and can provide better evidence for the presence 
of the HM in the waveforms. 



B. Extraction of Mass and Black Hole Spin 

As a second test, we will extract the mass AI and the 
spin parameter a from the waveform. We extract the 
mass first. To do this, we set a /Mo to the value given in 
[3^ 1 for the final black hole, but allow M to vary (note 
that the QNMs and the HM are function of a/M, not 
a/Mo). We then calculate Pmax as a function of M, us- 
ing equation (|66p , and find the value of M which maxi- 
mizes Pmax- When doing so, there is a distinct residual 
oscillation in the extracted value of M /Mo- Investigation 
reveals that the residual oscillation is compensated for by 
including the first £ = 4, m = 4 corotating QNM in the 
fit. Appendix C gives a brief discussion of the possible 
sources of this mode mixing in the numerical waveform. 

The top panel of Figure [5] gives the extracted M as a 
function of to, for two sets of frequency modes. The first 
set is composed of the first three ^ = 2, m = 2 QNM pairs 
(co and counterrotating), plus the first corotating i = A, 
TO = 4 QNM. The second set replaces the third £ = 2, 
TO = 2 QNM pair with the HM. We expect the extraction 
to fail at early values of to , where the sinusoids are a poor 
fit to the waveform, and at late values of to, where the 
waveform has decayed significantly. Indeed, one can see 
in the top panel of Figure [5] that the extraction begins to 
diverge as portions of the waveform preceding the peak of 
I'i/'l at tinax = 3953. 8M0 (shown with a solid vertical line) 
are included in the extraction, corresponding to values of 
to earlier than the peak. The middle and bottom pan- 
els of Figure [S] show the results of the same extraction, 
using the first four and five £ ~ 2, m ~ 2 QNM pairs, 
respectively, and comparing to extractions which replace 
the QNM pair with the largest decay with the HM. We 
see in these figures that the extraction can be carried out 
to even earlier values of to than the peak of ■0, but since 
these extractions include portions of the waveform which 
do not correspond to ringdown, we do not expect these 
early time extractions to be accurate. 

In all cases. Figure [5] shows that the substitution of the 
HM does not improve the extraction over the next most- 
damped QNM. We find the mean and RMS deviation of 
the extraction over an interval to /Mo = [3954,4074] for 
all three extractions. This interval covers a region of to 
that begins just outside the peak and continues until the 
extractions begin to diverge rapidly. We compare to the 
mass given in |32l | giving a difference in extracted mass of 
AMqnm = Mf- Mqnm for the fits that use only QNMs 
and AMhm = Mf- Mhm for those that include the HM. 
We present the values of AAf in Table I. We see that in 
all cases, the HM gives a comparable extraction. For the 
extractions that use a larger number of QNMs, we see 
that the RMS deviation grows. This appears to be due 
to the fact that the extractions with a larger number of 
modes diverge slightly earlier than those with n = 3 and 
n = 2 QNMs shown in the top panel of Figure [5) 

We carry out the same test for the spin of the black hole 
a /Mo, by fixing the mass M at its final value and allow- 
ing a/ Mo to vary. The results are given in Figure [51 The 



12 



n QNMs 


AMqnm 


AflQNM 




(Mo X 10"^) 


[Mo X 10-=^) 


3 


1.23 ± 1.66 


0.20 ±3.36 


4 


1.70 ± 1.96 


1.23 ±4.95 


5 


0.84 ±8.14 


2.55 ±8.74 


n QNMs 


AMhm 


Aohm 




(Mo X 10"^) 


(Mo X 10"^) 


2 


0.53 ± 1.37 


0.42 ± 2.95 


3 


1.15 ± 1.73 


1.57 ±5.07 


4 


0.40 ± 8.03 


3.00 ± 8.93 



TABLE II: Extracted masses and spin parameters, for extrac- 
tions using the first n QNMs, and extractions using tlie first 
n QNMs plus tlie HM. 



situation is the same as in the case of the mass extrac- 
tions. We calculate the mean and RMS deviation from 
the mean on the interval to/Mo = [3954,4074], and com- 
pare the extracted spin parameters to that given in [s^ . 
and give the results for Aa in Table IIVBI Once more, 
extractions with just QNMs are essentially the same as 
those with the HM replacing the most rapidly decaying 
QNM of a given set. From these extraction tests, we 
cannot conclude that the HM is present in the numerical 
waveforms. 



V. DISCUSSION 

Using two methods, we have found an additional ring- 
down mode for the Kerr black hole. This HM depends 
only on the fundamental properties of the black hole: 
it oscillates at the horizon frequency of the black hole, 
and decays at a rate proportional to the surface gravity 
of the black hole. It will arise when generic initial per- 
turbations source linear gravitational radiation, a situa- 
tion that would occur as the spacetime transitions from 
a regime of stronger, nonlinear perturbations into a fi- 
nal ringdown phase. This occurs at the last stage of a 
compact binary merger, or stellar core collapse resulting 
in a black hole. We emphasize that this mode is not in 
the QNM spectrum which is generally taken as the com- 
plete spectrum for the ringdown of a black hole. At the 
same time, this oscillation mode is part of what is nor- 
mally considered the "ringdown" phase of an event that 
results in a final black hole, since it arises in linearized 
perturbation theory about the final black hole. 

In fact we have discussed two possible decay rates for 
the HM, each dependent on our model of how the space- 
time transitions from nonlinear evolution into the linear 
regime. One mode, found using a naive model of tran- 
sition at a set time slice, decays rapidly. The second 
was found by noting that nonperturbative regions of the 
spacetime should be bounded by in-going and out-going 
characteristics, and is physically better motivated. It 
has a decay rate 7 = 2gH, approximately the same de- 
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FIG. 5: (color online) Values for the extracted masses Mqnm 
(solid) and Mhm (dashed) as functions of to, for three extrac- 
tions. The top panel compares the first three QNM pairs with 
the first two and the HM. The middle panel compares the first 
four QNM pairs with the first three and the HM. The bottom 
panel compares the first five QNM pairs with the first four 
and the HM. The solid horizontal line is at M/Mq = 0.95162, 
the mass of the final black hole as given in [s^]- The solid 
vertical line gives the peak of \^f}\. 
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FIG. 6: (color online) Values for the extracted spins sqnm 
and flHM as functions of to, for three extractions. The top 
panel compares the first three QNM pairs with the first two 
and the HM. The middle panel compares the first four QNM 
pairs with the first three and the HM. The bottom panel 
compares the first five QNM pairs with the first four and the 
HM. The solid horizontal line is at a/Mp = 0.65325, the spin 
of the final black hole as given in [32]. The solid vertical line 
gives the peak of \^\. 



cay rate as the n = 3, £ = 2, m = 2 QNM. We find that 
this mode has the same influence on overlap calculations 
as the n = 4 pair of QNMs, though it does not appear 
to improve parameter extraction over the use of an ad- 
ditional QNM pair with n > 4. Due to its comparable 
decay rate to the n = 3 QNMs, it should be considered 
in the construction of waveform templates that use n > 3 
QNMs. The HM should also be included as part of the 
ringdown spectrum when considering the potential use 
of an observed ringdown signal as a test of the No-Hair 
Theorem. Otherwise, the presence of this non-QNM os- 
cillation in the spectrum might lead one to conclude that 
the signal was emitted from an object other than a Kerr 
black hole. 

The analytic approach presented here also builds some 
intuition into the origin of various frequency modes of lin- 
ear perturbations of the Kerr spacetime. The HM studied 
here arises when the influence of perturbations near the 
horizon are considered. Integration of these initial per- 
turbations using the Green's function approach results in 
the presence of a pole in the frequency integral of equa- 
tion (|11|) . This mode depends only on the properties of 
the black hole which govern its near horizon geometry. 
Meanwhile, the usual QNMs arise because of the poles in 
the Wronskian of the radial Green's function, equations 
dm -([22]). In our model, these modes arise due to the 
interaction of the initial perturbations with the compli- 
cated potential of the wave equation present further from 
the event horizon, a situation analogous to that explored 
for Schwarzschild black holes by Price pGl ) . In this work, 
decaying perturbations on the surface of a collapsing star 
are associated with out-going radiation; comparison of 
our results with [26] indicates that our HM is associated 
with the decaying mode at the stellar surface, but that 
the rotation of the Kerr black hole in our case guarantees 
that this mode oscillates with the horizon frequency in 
addition to its simple decay. 

We have also reviewed the problem of gravitational 
radiation from a point particle infalling near the hori- 
zon. Previous work [l^ both motivated this study and 
guided our investigation. However, our results conflict 
with those of the motivating study. In investigating this 
discrepancy, we have found an error in the original cal- 
culation of |16i], the correction of which cancels the first 
order results for the radiation at infinity. We have also 
argued that the form of the next order correction agrees 
with our results for vacuum perturbations. We leave the 
detailed calculation of the correct second order terms to 
a future study jsTj . 

Future study using a variety of numerical waveforms 
will be key in determining the importance of the HM 
in template generation and gravitational wave detection. 
In a simulation where the excitation of slowly decaying 
QNMs is suppressed, we would expect the HM to be a 
clear component of the ringdown. Future study of how 
one might suppress this QNM excitation would be valu- 
able, and such simulations would provide the best testing 
ground for the presence of the HM in numerical Simula- 
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tions. In addition, the properties of the near horizon 
region, the HM itself, and the regularity conditions on 
the initial data discussed here may be of interest in the 
mathematical study of the stability of the Kerr black hole 
(see e.g. 



34| and the references therein). 
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Appendix A: Green's Function Formalism For the 
Teukolsky Equation 

The Teukolsky equation, for spin s = — 2, can be writ- 
ten in Boyer-Lindquist coordinates using the Kinnserly 
tetrad as [T^ 



(Al) 



with L the linear Teukolsky operator described below; 
■0 = '^i/p^; and the source term T a complicated 
function of the stress-energy tensor T**", the Kinnersley 
tetrad, and the Kerr rotation coefficients. The Teukolsky 
operator is 



L = Lr + Lg + Aidl + A^dt + A^dtd^ + Aidl 
+Azd^ + Ac , 

1 



Le = - 
Al = 

A2 = 
^3 = 
^4 = 

A, = 



90 (sin Ode) 

.2 „• 2 c 



A2 sin 61 

(r^ + a?Y " "^^^ 

A3 A2 ' 

4M(r2-a2) 4(r + wcos^ 

A3 
AMar 



(A2) 
(A3) 

(A4) 



A2 



A3 
2 



1 



A3 A2 sin^ e ' 
4a(r-M) icos0 



A. 



A3 A2 sin^ e ' 

4 cot^ e + 2 



A2 



(A5) 
(A6) 
(A7) 
(AS) 
(A9) 
(AlO) 



We introduce the adjoint operator L* , which is 
the Teukolsky operator with the substitutions {dt — > 



-d(f,), and the Green's function G{x''^;x'^) 



for L, which obeys L[G{x''']x'')] = L*[G{x'i']x 
6{t' - t)S{r' - r)S{e' - e)S{(t)' -4)) = S^ix'f' - x''). Now, 
note that given a pair of functions u and v we have 

uL*[v\ — vL[u] = dr[/^^^{vdrU — udrv)] 

+ — — -dg\v sin Odgu — usinOdgv] 

+Aidt[udtv - vdfu] - A2dt[uv\ 

+A^[dt{ud^v) - d^{vdtu)] + Aid^[ud^v - vd^u] 

-A^d^luv] . (All) 

Now, we let u = ip^x^), v = G{x'^; x^), and we in- 
tegrate the entire expression over the domain of interest 
for our situation, t £ [0,oo), r e (r+,oo), 9 G [0, tt], 
(j) G [0, 27r], using spherical polar coordinates and a Eu- 
clidean volume element d^x = dtd^x — r"^ smO dtdrdOdcf). 
To evaluate the left side of equation (|A11[) we note 



d^x '^|J{x^')L*[G{x'^'■x^')] = J d'^x ^ix'')S\x''' ~ xf") 

= tPix'^") . (A12) 

Also, we have 

J d'^x G{x'^';x'')L[^p{x^')] ^ J d*x G{x'^']x^')/^-^T . 

(A13) 



On the right hand side of equation (jAll[) . we note 
that the terms involving ^4 and ^5 vanish when inte- 
grated over (j), due to the periodicity of (j). The term 
—A3d^{vdtu) vanishes for the same reason. The first 
term becomes a boundary term when integrated over r. 
In order to have only in-going waves at the horizon, and 
only out-going waves at infinity, we must impose homoge- 
neous boundary conditions on tp{x'^) and G{x'^;x'^), and 
so this term also vanishes. The term involving derivatives 
of 9 vanishes when integrated over 9, since we require that 
the initial data and the Green's function be regular on the 
boundary of [0,7r]. The terms involving Ai and A2 are 
total derivatives in time, and so when we integrate over 
t wc remove the time derivatives and evaluate the terms 
on the boundary at t — 0. Since our physical source is 
transient, the terms vanish at the bound of i — > cx). We 
have 



?/'(x'^) 



d'^x G{x''';x'')A-^T 
+ I d^x AidtGix'^';x^')ip{xf') 
d^x AiGix'^';x^')dt^p{x^') 
d^x A2Gix''';x'')'iljix^') 



t=0 



t=0 



d^x A3^^Gix'^'■,x^')'^lJix^') 



(A14) 



t=o 
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Thus far we have kept the source term T in place for 
comparison with other studies of the evolution of the 
Teukolsky equation. In this study we are interested in 
the vacuum case, and so we set T = here and through- 
out Section II. 

We can further simplify this expression for 7p{x^) by 
performing the angular integrations. Let us expand the 
initial perturbation in terms of spherical harmonics, 



t=0 



In addition, we expand the Green's function in the fre- 
quency domain, where it can be written down explicitly 
in terms of the spin-weighted spheroidal harmonics and 
the Green's function for the radial Teukolsky equation 



G(x'^;x'^) 



(A17) 

So, we have ^^G{x'^';x^') = -iTOG(x'^; x^). With this, 
we can now perform the integration over </>, using the 
identity 



27r 



27r 



i(m-m')0 _ r 



(A18) 



which allows us to resolve the summation over m' con- 
tained in equation (|A15P - (|A16[) . From here, it is con- 
venient to impose the near horizon approximation, for 
which the motivation is discussed in Section IIB. We keep 
terms only to the leading order in e = (r — r+)/r+ <^ 1. 
In this approximation, we have that A 2A/r+Ke, with 
K = -^/l — a? /A'P. To first order in e. 



-imA-i 



Ax 

A2 



{2Mr+)-^{Ke)-^ 

2Mk + ima 
~2(Mr+)2(Ke)3 ' 



(A19) 
(A20) 



We note that all 9 dependence for these functions enters 
in at second order in the near horizon expansion. We 
define 



(2£' + l)(f -m)! 
4Tr{£' + my. 



X / sm9d0Ptm{cose)SemUd)XA21) 
Jo 



(2£' + l)(r-m)! 
47r(f -l-m)! 



X / siu9dePtm{cose)Sim^{e),{k22) 
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where Pim{x) are the associated Legendre polynomials. 
The functions a(jn^{r) and Pimujif) are nonzero only on 
the interval r G [r+, (1 -I- ^)r4-], which allows us to trun- 
cate the radial integrals in equation (IA14p . In fact, we 
only desire the leading order behavior in e of these func- 
tions, and this is discussed in Section IIB. 

Inserting equations (|A15I) - (|A21I) into (|A14I) . and ex- 
changing primed and unprimed labels, we have finally 



/ 1^ E e"'"*+™^i?£™^(0^^™c.(e), (A23) 



dr' 



2Mr+{Ke) 



{2Mk + ima)aernui{r') 
2(Mr+)2(Ke)3 



Gemujir,r') . 
(A24) 



We resolve this expression in Section IIB. 



Appendix B: The Teukolsky Equation in the 
Newman-Penrose Formalism 

We refer the reader to [13, HI] for the full formalism. 
Here we simply collect some of the longer expressions 
used for Section IIIA. 

From (2.14) of [l3| we have 



(A 37 - 7 4^ -I- p.){D + 4e - p) 
-{S - f + P + 3a + Att){S - t + 4/3) 



-3*2 



*f = 47rr4. (Bl) 



Here refers to the background value of the NP scalar, 
^'2 — Mp^ for Kerr. The scalar is the perturbative 
value of ^'4, which is zero at leading order for Kerr. Here, 
D, A, S are all derivative operators along the directions 
of the null basis, and the Greek characters represent com- 
binations of the spin coefficients. Also note the unfortu- 
nate but standard use of tt on the left-hand side to refer 
to one of the spin coefficients in the null tetrad, while on 
the right side it refers to the numerical tt from the Ein- 
stein field equations. It is generally clear which is which, 
and in any case the NP coefficient enters at subleading 
order here. The source term T4 is given by 

T4 = (A + 37 - 7 + 4^^ + /2) [(J - 2f + 2a)r„^ 

-(A 27- 27 + /i)T™ 

+ {S - f + ^ + 3a + Att)\{A + 2"f + 2^)r„,-n 

-{5-f + 2$ + 2a)Tnn\ .(B2) 

Here, the terms Tat are the components of the stress- 
energy tensor in the tetrad basis, r„„ = T^j^^n^n'^ , Tnm = 
Tf.^nf'rh'', etc. 
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To specialize to the near horizon approximation, we 
note that D contains A~^, and therefore dominates over 
all the other terms. In addition, we have 7 = 7 = pp{r — 
M)/2 and ^ = 0, to first order. 

Using the commutation relation between D and A 
(NP4.4), we have, near the horizon 



AD — D/S. — 2jD + (lower order terms) 



(B3) 



which subsequently gives the 0(A ^) term on the left- 
hand side of the equation 

{DA + 47Z))*4 (B4) 

We investigate this expression more fully in Section IIIA. 

Appendix C: Mode Corrections to the Wavefunction 

The presence of oscillation at the n=l, i = 4, rn = 4 
corotating QNM merits some brief discussion. Figures [7] 
and [8] give the extraction of M and a using the first three 
QNM pairs (and comparing the the first two QNM pairs 
with the HM), without the £ = 4, m = 4 mode included. 
Comparison with the topmost panels of Figures [5] and [S] 
shows that the distinct oscillation is successfully removed 
by including this mode. 




0.948 



3940 3960 3980 4000 4020 4040 4060 4080 

to/Mo 

FIG. 7: (color online) Extraction of the mass M/Mo as a 
function of to, using only first three QNM pairs (solid) or the 
first two QNM pairs and the HM (dashed). Here we do not 
include the corotating n — 1, £ — 4, m = 4 mode. 

A certain amount of mode mixing between the QNMs 
is expected due to the fact that the waveform is decom- 
posed into spin- weighted spherical harmonics during the 
extraction of the waveform. In fact, the angular eigen- 
functions of the Teukolsky equation are the spin- weighted 
spheroidal harmonics. These functions become the usual 
spherical harmonics when aw = 0. Using this fact, the 
spheroidal harmonics can be expanded in terms of spin- 
weighted spherical harmonics and powers of aco, as first 



discussed in [2C| . Only spherical harmonics with the same 
s and m contribute in the expansion. As such, we see 
immediately that the mixing with the £ — 4, m — 4 
0.660 [ 
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FIG. 8: (color online) Extraction of the spin a/Mo as a func- 
tion of to, using only first three pairs QNMs (sofid) or the first 
two QNM pairss and the HM (dashed). We do not include 
the corotating n = 1, 1 = 4, m = 4 mode. 



QNM frequency cannot arise from the decomposition into 
spherical harmonics. The portions of the waveforms that 
can mix into the ^ = 2, to = 2 waveform arise from the 
expansions of 5'32,S'42, etc. Explicitly, the expansion of 
the spheroidal harmonic for s — —2, is 



S^m — 



2Yim 



4auj 



21' + 1 + 1) 



+0(a2^2 



f (^' + 1)] 
(CI) 



where C^^^^ are the usual Clebsch- Gordon coefficients. 
For both the horizon mode and the lowest order QNMs, 
ao; < 1 is true for all a/M, and so the expansion is 
not obviously divergent, although it is only good when 
a/M ^ 1. The inclusion of additional QNM frequencies 
with TO = 2 does not remove the residual oscillation in 
the extraction of M and a seen in Figures[7]and[5] (though 
a corotating ^ = 3, to = 2 reduces the amplitude of the 
oscillation somewhat). In fact, extractions using a large 
number of modes generally have sharp features, in addi- 
tion to systematic deviations from the values of M and 
a given in [32 1. 

The presence of the I = 4, m = 4 mode in the 
£ = 2, m = 2 waveform is unexpected, and we attribute 
it to errors arising from the numerical generation and ex- 
traction of the waveform. The spectral code used in (32j 
generates its gauge dynamically, and while the waveform 
extraction method attempts remove gauge effects, studies 
find that these gauge effects still generate errors [35| . We 
suspect such gauge errors are the source of mode-mode 
mixing. 
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